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!•  Introduction 

^r!5ie  major  outstanding  probleirsin  the  theory  of  the  pulsating 
gas  bubble  produced  by  an  underwater  explosion  is  the  deviation 
of  the  bubble  from  the  spherical  shape,  and  the  effect  of  this 
deviation  on  the  secondary  pressure  pulse*  As  a  step  to  the 
solution  of  this  problem,  we  shall  suppose  that  the  gas  bubble 
Is  constrained  to  be  ellipsoidal  in  shape.   The  flow  of  the 
water  due  to  a  pulsating,  moving  ellipsoid  can  be  determined  ex- 
plicitly in  terms  of  elementary  functions  and  leads  to  a  system 
of  differential  equations  for  the  behavior  of  the  bubble.   These 
differential  equations  are  derived  in  this  note,  and  they  show'^he 
deviation  from  spherical  shape  is  very  marked  when  the  bubble  is 
in  its  contracted  stage  near  the  minimum  size.   More  exact  informa- 
tion awaits  the  numerical  solution  of  these  equations , which  is 
under  way. 

The  reason  for  a  change  in  shape  from  the  spherical  is  the 
following.   The  bubble  in  expanding  and  contracting  imparts  a 
linear  momentum  to  the  surrounding  water,  due  to  the  buoyant  action 
of  gravity  and  the  presence  of  nearby  surfaces.   This  causes  the 
bubble  to  move,  and  because  of  conservation  of  momentum,  the  motion 
is  rapid  when  the  bubble  is  near  its  minimum  size. 

It  is  a  result  of  classical 
hydrodynamics  that  a  moving  sphere 
has  a  high  pressure  at  the  nose 
and  rear  of  the  sphere  and  a  low 
pressure  along  the  sides  (see  the 
diagram).   Thus  the  gas  bubble  in 
moving  will  be  flattened  out  in 
the  direction  of  motion.   This 
flattening  process  can  be  followed  by  assuming  the  bubble  to  be 
ellipsoidal.   (Of  course,  it  is  quite  possible  for  the  change  in 
shape  to  become  more  drastic  than  the  ellipsoidal;  e.g.  the  shape 
might  become  rather  wavy  or  it  might  actually  become  toroidal. ) 
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II.  The  kinetic  energy  of/ water. 

1.  To  describe  t^e  ellipsoidal  bubble,  we  shall  use 

spheroidal  coordinates.   See 
Lamb,  Hydrodynamics,  6th 
edition.  Chap.  V,  esp.  p. 142-6. 
The  bubble  is  assumed  to  be  an 
oblate  spheroid  (planetary 
ellipsoid),  with  semi -axes 

c^  ,  «i  ,  ^  where  ^  <  oc  • 
The  Intersection  of  the  zx-plane 
with  the  spheroid  is  an  ellipse 
with  foci  in  the  x-axis  at  a  distance  kf  rom  the  origin,  whefe 

r~o     ^ 

k  =  y oc   -  v>      •   See  the  diagram.    Spheroidal  coordinates   are  introduced 
in   space.      These  are  coordinates    (  ^  y  lA  ,  (S>  )   vh  ere  0  -^  -'^  t 
-1     5yU.<     1,    0SG^5     2Tr   and  where     x,y,  z  are  related  to  ^  ,//<.,  <25 
through  the  equations: 


y  1  +^  ^  yi  -  /x^cos  ^  0  5^< 


x  =  kYl+4Vl-  ^    cos  CQ  0  5  4    <   c>o 

2./,        ..Z. 


(1)  <y  =  k  ]j  1  +5"     y  1  -A     slnCp  0  <  i/  ^     1 


z  =  k  ^u.  0  ^  (^  ^     2iT 


where  k  is  the  semi -focal  distance  of  the  spheroidal  bubble. 
The  coordinates  r>     ^  /^   are  more  eustomarllY  presented  da  ^=  sinh^ 
jX-   cos  ©  J  then^l  +^   =  coshvi  ,|/l  -^x   =  sin  9   •   The  coordinate 

surfaces  Z  =  constant  are  a  family  of  confocalel lips olds,  with 

/ 2*  '  2         -y 

semi -focal  distance   k  and  semi -axes  k     ^  \  -^  ^    ,k4/l+Z    ,    k^ 

respectively. 

In  particular,    the   equation  of  the  surface   of  the    spheroidal 
bubble   is 


t- 


^    ^ 


^         (a   constant). 
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Thus,  the  spheroidal  bubble  is  characterized  by  two  quantities 
namely  k  and  <^  •  (The  semi -axes  o^    9    '<  ,    y    will  no  longer  be 
used.)   The  coordinates  (^  ' /^  *   ^^  defined  by  (1)  are  then 
the  spheroidal  coordinates  of  an  arbitrary  point  in  space. 
The  element  of  length  dS  in  space  expressed  in  shperoidal 
coordinates  is 

(2)  ^2  ^  l£(£fi^  ^2  ^  k^(£!±4i  dz.2  ,  k2(i  ^  ^2^(1  .^2) 

2.   We  shall  now  consider  the  flow  produced  in  the  water 
surrounding  the  ellipsoidal  bubble  when  the  ellipsoidal  changes 
Its  shape.   The  change  in  the  ellipsoid  can  be  produced  by 
changing  either  k  or  4  «»  ^''^^  ™®  shall  first  suppose  that  k 
changes,  while  ^   is  fixed.   This  correspondlaae  to  a  similarity 
transformation,  and  is  a  "pure"  expansion.   If  we  follow  a 
point  on  the  surface  of  the  ellipsoidal  bubble  with  given 
coordinates  (^q>/^>  ^  ^»  ^^®^  *^®  velocity  of  that  point  Is 
by  (1):         ^ 


2 


(3)  ^  y  =  kVr77^1/rr;?sinc^ 

5  =  ^  ?'  o/^ 

On  the  other  hand,  the  direction  of  the  normal  to  the  ellipsoid 
at  the  point  (^o>/*»*^)  is  given,  also  from  (1),  by 


r, 


9i  =  k-«L2^  V^ 


a7  ^  ^>,7r*T*  VI  -AC  coso?  /direction 
^    Vl  +<„     /      ^   normal  to 


(4^  \  ellipsoid  at  . 
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The  component  V   of  the  velocity  of  the  surface  In  the 
direction  normal  to  the  ellipsoid  is 


(5)  Y.  =  "" 


ImUJlU:!! 


'V#^^S^ 


A  simple  ealculEtlon  using  {3)  and  (4)  gives 


The  boundary  condition  for  the  potential  function  $^ 
describing  the  flow  is 

\1)  -^  =  V   on  5^=  t^ 

3n     n     -'    ^  o 


By  (2), 


~2 


(8)  -—i  =  (-TT^)    or  <  =  ? 


v^-^ 


Combining  (6),  (7),  (8),  we  therefore  have 
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Thus, y^     is  constant  on  the  surface  of  the  ellipsoid, 

and  the  potential  function^.  Is  therefore  easy  to  construct. 
We  have  in  fact 

where  f?-  (^)  =  cot"  ^^  ,  and  the  prime  on  q^  C^^)  means  the 
derivative  with  respect  to  ^ ,  (see  Lamb,  Chap.  V,esp.  page  143 
for  rotationally  symmetric  spheroidal  harmonicsji 
These  harmonics  have  the  form 

f.y-)  On'?' 

where  Pj-^/*"^  ^^^   Legendre  polynomials,  and  qj,(^)  are,  except  for 
multiplying  factors  making  them  real,  Legendre  functions  of  the 
second  kind  for  imaginary  arguments.   The  first  few  of  these 
functions  are 

P^U)  =  1  q„(^)  =  cot-^^ 


r 


(10) 


-1. 


P^^)  =LL,  q^(^)  =  1  -/cot'-^^; 


3.   Consider  the  change  in  the  ellipsoidal  bubble  due  to 
varying  ^  ,  keeping  k  fixed.   Following  a  point  on  the  surface 
of  the  ellipsoidal  bubble  with  f ixed  t*.  ,  QP  ,  we  have  for  the 
velocity  of  that  point,  by  (1): 


ro?o 


:1^ 


X  =  k  ' '  '  '   ■ 1/1  -  ul.   cos  (Z? 
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The  component  ^„  of  the  velocity  of  the  surface  in  the 


n 


direction  normal  to  the  ellipsoid  is,  by  (4)  and  (5), 


^n  =  ^o"^) 


The  boundary  condition  for  the  potential  function  ^g  describing 
the  flow  is,  by  (7)  and  (8) 


-2    2 
2?o  *A 


ij'    I 


or 


?=  $ 


This   condition  can  be  written  in  the  following  form: 


a^   ■   ,.2 
b4       ^^ 


1-? 


2 
3 


^  (  f  ^2)^)  + 


)  or  /=  ^, 


3^^  -  1 


where   PgO^)  =  -^^- — --"^  •  1"^®  solution  for  ^g  is 


^o^^ 


(11)  ^P=-^ 


2         2, 


.2 


-^2(^0) 


Klo^^o) 


4.  Finally,  consider  the  ellipsoidal  bubble  as  fixed  in  size 
and  shape,  but  moving  in  the  z-direction  with  velocity  B  •  (The 
quantity  B  represents  the  vertical  distance  of  the  center  of  the 
bubble  from  same  level* )  For  the  velocity  of  a  point  on  the  surface 
of  the  bubble  moving  vertically  upward,  we  have 

^2  =  0,  .^=0,   ^  =  B  . 

at   'at     at 
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The  normal  coinponent  of  the  velocity  Is,  by  (4)  and  (5), 


By  (7)  and  (8),  the  boundary  condition  for  the  potential  function 
The  solution  for  ^^   Is 


(12)  #3  =  ^  ^    /x,qi(^)  . 


5.  The  potential  function  <±>  for  the  combined  flow  of  the 
water,  due  to  t he  changing  shape  of  the  ellipsoid  and  its  motion 
throughitthe  water,  is  the  sum  of  the  individual  parts  ^_,^g,^_ 
given  in  (9),  (IL),  (12).   It  is  more  simply  expressed  if  the 
volume  V  of  the  bubble  is  introduced.   The  volume  V  is 

The  rate  of  change  of  colume  V  is  simply  expressed  in  terms  of 
£:  and  ^  . 

The  potential  function  <^=  ^^  +  ^2  +  ^^   for  the  combined 
flow  is,  from  (9),  (11),  (12), 


(13)     f>=JL  -  .^.     B  k 
4Trl 


?k  '^o(^)  ^         ,     ^q  (^) 

-Q-,(r^)'    ^ 


Ql^^o' 


to^^^ 


3(l+ro^(-q;(<o)^ 
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The  convenient   form     -j— r  Q   (Z)   for  the    first   terra  Is   a  matter 

4trK:  ^o  ^ 

for  simple  calculation,  but  It  can  also  be  obtained  by  the 
follov/lng  general  argument.   At  a  large  distance  R  from  the 
center  of  the  bubble,  the  potential  function  must  by  conservation 
of  mass  be  of  the  form 

At  large  distances  from  the  origin 

1o'<'    1    1 
k      k    R  ' 

whereas  all  the  other  terms  In  (13)  vanish  to  a  higher  order  In  —  . 

R 

6»   The   kinetic   energy  ^  of  the  water  can  now  be   obtained. 
We  have  by  Green's   formula, 

T    =  I  P///(grad^)2dr    =  | ///$(-  •|i)dS   , 

where  j>     is  the  density  of  the  water  and  the  Integral  Is  extended 
over  the  surface  of  the  ellipsoid.  Using  (8)  and 


dS    =   kl^^' 


^sJZL  d^    .ky(l+^^)(l   y-2)dcy 
from  (2),   we  find 


(14) 

1 


=  Try'k(l  +^^)  /^(.^)d;u- 
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Prom  (13), we  have 


a|« 


at 


b        4TTk(l    +^^) 


+  BkJi^  + 


r 


5(1+^) 


Po(/^)     . 


2^ 


Substituting    (13)   and    (15)    into    (16)   and  using   the   orthogonality 
of  t  he  Legendre  polynomials   P^^Cw-),   we   find 


exr 


:v^ 


The  quantity  k  can  be  V  and  ^  .  Instead  of  V  we  shall  use, 
however, the  radius  A  of  the  sphere  with  the  same  volume.   Thus,  set 


(17) 


V=|.A^ 


Ther  elation  between  k  and  V  yields 


k  = 


A 


<y'(i*?!)V3 


These  will  be  substituted  in  (17).   But  one  final  simplification 
the  subscript  o  on  ^  will  henceforth  be  dropped.   The  result  is 


where 
(19) 


o 

3»2 


IaS^  2 


m^(^)  =  5^3(1  4^2jl/3  q^(^) 


m. 


1^      ^ 


=^(4)  = 
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-lo- 
in the  formulae,  the  size  and  shape  of  the  ellipsoidal 
bubble  are  described  by  the  following  two  quantities: 

(1^  The  radius  A  of  thesphere  with  the  same  volume  as  the  ellipsoid. 

(2)  The  quantity  ^  which  determines  the  sJiape  of  the  ellipsoid 
and  represents  the  ratio  of  the  semi -minor  axis  y  to t he  semi -focal 
distance  k, 

^    k 
In t  erms  of  the  semi -major  and  semi -minor  axes    and    , 


OC 


/w 


The  position  of  the  ellipsoid  is  determined  by  the  vertical  distance 
B  of  its  center  from  some  horizontal  level. 

The  quantity  ^    can  vary  from  0  to  «>o  ,  the  case  Z    =  »<=> 
corresponding  to  a  sphere  and  thee  ase  ^  =  0  to  a  flat  circular 
disc.   For  the  spherical  case,  Z^     =  ©o,  we  have   m  =  1,   m,  =  1, 
Kg  =  0  and  the  kinetic  energy  formula  (19)  reduces  to  the  known 
spherical  case. 

In  view  of  (18),  it  is  convenient  to  introduce  a  new  variable 
A  in  place  of  ^  .   Set 


(20)         \  =   /  /^^)dr 

As  X     varies  from  <>«  to  0,  X  varies  from  0  tothe  finite  value 

/  y^C'T^)  d^   t    the  value  A  =  0  corresponding  to  a  sphere. 
?En  t  erms  of  A  the  kinetic  energy  is 

(„>)     X  =  ,  .A^l^  ^      i  ^S|2  ^  1  ^ty2 
^  2m  ^  -L   6        2 

where  m^,  m^  are  functions  of  X  defined  through  (19)  and  (20). 
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III.  The  equations  of  motion  of  the  bubble* 

1.  The  potential  energy  U  of  the  gas  bubble  Is 
U  =  Pg(D  -3)  Itta^  +  G(A) 

where  D  as  the  Initial  depth  of  the  explosion  from  33  ft.  above 
sea  level,  B  Is  the  upward  displacement  of  the  bubble  at  any  time, 
and  G(A)  its  the  Internal  energy  of  the  gas.   The  terra  >o  g(D-B) 
represents  the  hydrostatic  pressure  at  the  center  of  the  bubble. 
The  energy  equation  is 

'y  +  U  =  E  (  a  constant)* 

Introducing  dlmensionless  variables  as  in  the  AMJ-NYU  report 
"The  best  location  of  a  mine  near  the  bottom  of  the  sea"  by 
Friedman  and  Shiffman,  the  energy  equation  becomes 


T  u 

In  addition  to  the  Lagrange   equation  for  b,    there   is   a  Lagrange 
equation  for    X  •      These  equations   are 


a.  (™  .  l,.3f,l     =  i-  a' 


O 

loA\                                 d      /    5;  ,        "^o        3«2  %  1     3^2 

(24)                   -          --—  (a    A)   =  -T^  .a  a  +  — ^  •  r  ^  ^ 

dt                          dX  dX  6 


The  equation  (22),  (23),  (24)  are  the  differential  equations  for 
determining  a,  b,  X  as  functions  of  t. 

2.  The  Interest  centers  around  the  contracting  stage,  especially 
when  the  bubble  is  near  Its  minimum  size.   At  this  time,  equation 
( 23 )  shows  that  the  linear  momentum  ' 

(25)  m^  .  -i  a^b  =  s 

1   3 
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Is  practically  constant  and  equal  to  s,   which  can  be  calculated 
in  terms  of  d  .   In  case  there  are  nearby  *igld  walls  or  free 
surfaces,  this  conclusion  is  still  valid,  except  that  s  now 
depends  on  the  distances  to  the  nearby  objects.  Substitution  In 
(22)  and  (24)  gives  approximately 

-2 

(26)  m  oaV  +  -H^     +  h^ '^^  +  a^  +  -^  =  1 

2m,a  a  / 

J        c  •           <im          -z   r>        1      dn,  „-2 

(or/\  d   /„5  -x  » o      „3.2    .    1     1  OS 

(27)  -(a    X)    =_   .a  a     +  ^ -[>      '  ^ 

These  equations  are  still  approximately  valid  throughout  the 
entire  contracting  state,  since  for  larger  a  the  errors  made  in 
using  the  equation  (25)  are  small* 

In  the  case  whete  there  are  nearby  surfaces,  the  only 
appreciable  effect  on  the  equations  (26), (27),  besides  the  influence 

on  the  value  of  the  momentum  constant  s,  is  an  extra  factor  (1+f  ) 
occ\irrlng  in  the  terms  a  S^  in  (26)  and  (27).   All  other  effects 
are  negligible. 

5.  An  analysis  of  equation  (27)  shows  the  extreme  instability 
of  the  shape  of  the  bubble.   Near  the  minimum  size,  the  "shape" 
momentum 

a 


5x=  A 


has  a  value  given  by  the  time  integral  of  the  right  hand  aide  of 
(27).   Then 

which  shows  how  large  /.   can  become  for  small  a.   Thus  X  may 
become  large  and  near  its  upper  limit,  which  means  a  flat  ellipsoid. 

In  this  case  m  is  small,  and  m,  is  large.   Also  both  g^  and 

I  dm, 
;;pi— Tn'   are  small.  Exact  statements  about  the  size  of  A  require  the 

numerical  integration  of  (26)  and  (27). 
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4.  The  displacement  of  the  bubble  can  be  determined 
from  (25)  by  solving  for  b 

•i   ss 


and  integrating  with  respect  to  t«   This  shows  that  for  large  m- , 

'^   b,  Ab  and  the  translational  energy of  the  water  are 

T  -.  2m-  a"^ 

small.  1 

The  pressure  pulse  P  emitted  by  the  contracting  bubble  is 
determined  from  Bernoulli's  equation 

where  the  remaining  terras  are  dropped  because  they  aa©  down  at 


a  faster  rate  than  1/R.  Using  the  first  term  of  (15),  which  is 
dominant  for  large  distances,  we  have  in  dimensionless  form 

P  =  (a^a)- 

Differentiating  (26)  with  respect  to  t,  we  see  that  the  peak  pressure 
p  is  expressed  with  the  coefficient  m  in  the  denominator.   Thus, 
the  smaller  m  is,  the  greater  as  the  peak  pressure  p  of  the 
secondary  pulse. 

All  of  the  remarks  in  th«S  section  apply  to  any  shape  bubble. 
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